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Abstract We present a novel algorithm for calculating the discrete fractional 
Hadamard transform for data vectors whose size A is a power of two. A direct 
method for calculation of the discrete fractional Hadamard transform requires 
multiplications, while in proposed algorithm the number of real multipli¬ 
cations is reduced to Alog 2 A. 
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1 Introduction 

Discrete fractional transforms are the generalizations of the ordinary discrete 
transforms with one additional fractional parameter. In the past decades, var¬ 
ious discrete fractional transforms including discrete Fourier transform 
discrete fractional Hartley transform [3] , discrete fractional cosine transforms 
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and discrete sine transform |4] have been introduced and found wide applica¬ 
tions in many scientific and technological areas including digital signal process¬ 
ing 0, image encryption 0, 0, 0 and digital watermarking [S] and others. 
Different fast algorithms for their implementations have been separately de¬ 
veloped to minimize computational complexity and implementation costs. In 
[To] a discrete fractional Hadamard transform for the vector of length N = 2^ 
was introduced, however a fast algorithm for the realization of this transform 
has not been proposed. In our previous paper m we describe a rationalized 
algorithm for DFRHT possessing a reduced number of multiplications and 
additions. Analysis of the mentioned algorithm shows that not all of existing 
improvement possibilities have been realized. In this paper, we proposed a 
novel algorithm for the discrete fractional Hadamard transform that require 
fewer total real additions and multiplications than our previously published 
solution. 


2 Mathematical background 


A Hadamard matrix of order N is a N x N symmetric matrix whose entries 
are either I or —I and whose rows are mutually orthogonal. In this paper we 
will use the normalized form of this matrix and we will denote it by Hjv. For 
A^ = 2" the Hadamard matrices can be recursively obtained due to Sylvester’s 
construction m- 


H2 


I 

71 


I 1 
I -1 ’ 


H^ = 


1 

71 


Hn h« 

2 2 

Hn -Hiv 
2 2 


( 1 ) 


for Af = 4,8,...,2". 

Definition of the discrete fractional Hadamard (DFRHT) transform is based 
on an eigenvalue decomposition of the DHT matrix. Any real symmetric ma¬ 
trix (including the Hadamard matrix) can be diagonalized, e.g. written as a 
product m 


N-l 

= Z^AjvZ^ = ^ Afcz7(z7)^ (2) 

fc =0 

where A^v is a diagonal matrix of order 2", whose diagonal entries are the 
eigenvalues of H^v 


Ajv = 


Aq 

•^1 

0 


0 


Aat-i 


( 3 ) 


Zat = [z7 I z7 I ■ • ■ I ^7 matrix whose columns are normalized 
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(k) 

mutually orthogonal eigenvectors of the matrix Hv- The eigenvector is 
related to the eigenvalue . A superscript T denotes the matrix transposition. 

The DFRHT matrix of order = 2" with real parameter a was first 
defined in m- This matrix can be regarded as a power of the DHT matrix, 
where the exponent a = al'K 


N-l 


H a _ 

N — 




fc =0 


(fc)/ (k)-^T 

z- (ziv ) 




(4) 


For a = 0 the DFRHT matrix is converted into the identity matrix, and for 
a=l it is transformed into the ordinary DHT matrix. Generally the DFRHT 
matrix is complex-valued. 

An essential operation, by obtaining the discrete fractional Hadamard ma¬ 
trix, defined by Q, is calculating the eigenvalues and the eigenvectors of the 
matrix Hjv. The only eigenvalues of the unnormalized Hadamard matrix of 
order A = 2" are known to be 2"/^ and —2”/^ [ 13 ], hence the normalized 
Hadamard matrix H^v has only the eigenvalues 1 and —1. A method for find¬ 
ing the eigenvectors of Hadamard matrix was firstly presented in m, but in 
[10| a recursive method for calculation the eigenvectors of the Hadamard ma¬ 
trix order 2”+^ based on the eigenvectors of the Hadamard matrix of order 
2"’ has been proposed. We will use this method to obtain the DFRHT matrix. 
Here we will present it briefly. 

In [15] it was proven that if (k = 0,l,..., N — 1) is an eigenvector of 
Hadamard matrix of order A = 2" associated with an eigenvalue A, then vector 


-'(fe) 

^2N = 


(fc) 

{V2-l)vP 


(5) 


will be an eigenvector of the matrix H 2 Ar associated with the eigenvalue A. 

In [TU] it was proven that if is an eigenvector of H^v associated with an 
eigenvalue A, then the vector 


''2N ~ 


(1 - 72)' 

^(fc) 


AD 


( 6 ) 


will be an eigenvector of the matrix H 2 Ar associated with the eigenvalue —A. 
These two results allow as to generate the eigenvectors of Hadamard matrix of 
order 2"+^ from the eigenvectors of Hadamard Matrix of order 2". Knowing 
the straightforward calculated eigenvectors of the matrix H 2 


V 


(0) 

2 


1 

V^- 1 




1 - V2' 
1 


(7) 


associated with eigenvalues 1 and —1 respectively, the eigenvectors for Hadamard 
matrix of arbitrary order A = 2" can be recursively computed. In uni it was 
also shown so this recursively computed eigenvectors of matrix Hjv will be 
orthogonal. It should be noted that for any N = 2'^ there are only two distinct 
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eigenvalues of Hadamard matrix, so for > 4 the eigenvalues are degener¬ 
ated. Because of this fact the set of eigenvectors proposed in [15] and [T0| is 

(k) 

not unique. The igenvectors for k = 0,1,...,1V — 1, which are columns 
of the matrix Zjv (after normalization), as well as their associated eigenval¬ 
ues Afe, can be however ordered in different ways. In nni it has been also 
established a method of ordering the eigenvectors. In many cases, including 
the case of discrete fractional transforms is used so-called sequency ordering of 
the eigenvectors. This means that the fc-th eigenvector has k sign-changes. The 
discrete Hermite-Gaussians, eigenvectors of discrete Fourier transform matrix 
are ordered this way as well [2|. We will show this method of ordering of the 
eigenvectors in Example]^ 

Example 1 The number of sign-changes in eigenvectors and of matrix 
H 2 , determined by Q, is equal to 0 and 1 respectively. Using expressions © 
and Q we obtain the eigenvectors of matrix H 4 : 


■ 1 ■ 


■ -b' 


■ -b ' 


'9' 

b 

^( 0 ) _ 


v(i) - 

1 

^( 1 ) _ 

-b 

_ 1 

^4 — 

1 

b 

^^4 — 

1 

(M 

_ 1 

^^4 — 

-b 

1 


where b = \/2— 1. The numbers of sign-changes in the above vectors are 0, 1, 
3, 2 respectively (5 > 0). Therefore, a sequency ordered set of eigenvectors of 
matrix H 4 will be as follows: 


^(0) _ ^(0) _ ^(0) ^(2) _ -(1) (3) _ - (1) 

V4 — V4 V4 — V4 V4 — V4 V4 — V4 . 


The corresponding eigenvalues will be equal to: 


Aq — 1 Ai — —1 A2 — 1 A3 — —1. 

The relations obtained in Example can be easily generalized as follows: 


for 1 = 0 , 1 ,..., ^ — 1 and 


.(40 _ 

v(20 

2N — 

^2N 

,(4/+l) 


2N 

~ ^2N 

,(4/-S2) 

_ 4.(2i + l) 

2N 

~ ^2N 

,(4/+3) 

^(21 + 1) 

2N 

— ^2N 

Afe = 

{-1)\ 


( 8 ) 


(9) 


for fc = 0, l,...,2iV- 1. 

Both the eigenvectors of the matrix H 2 and the eigenvectors obtained for 

II (k) II 

higher order Hadamard matrices are not normalized. Let the notation 













Novel DFRHT algorithm 


5 


_ __ 

means the Euclidean norm of vector v]y . In [IT] it was shown that for any 
N = 2"^ we have the relationship 


dfc )||2 

'n II 


= (1 + 


( 10 ) 


where A: = 0, 1 , ..., iV — 1 and b — \f2 — 1. If we take the designation c= 1 + 6^, 
then normalized eigenvectors of Hadamard matrix of order N = 2'^ will take 
the form 


z 


(fc) 

N 





( 11 ) 


Using the normalized and sequency ordered eigenvectors of the Hadamard ma¬ 
trix, the eigenvalue decomposition ([^ of the Hadamard matrix can be written 
as follows: 



Hat = ZatAatZ^ = —VnAnVJj 
c 

(12) 

where Apj is 

the diagonal matrix whose non-zero elements are 



Afe = (-1)'= = 

(13) 

for fc = 0,1,, 
the form: 

...,7V— 1. Hence the definition Q of DFRHT matrix will take 

= \y^A%YI (14) 

where 

\a _ -jkTra 

(15) 


for fc = 0,1,..., iV — 1. 

Our goal is to calculate the discrete fractional Hadamard transform for an 
input signal xn in which the number of samples is equal to A^ = 2". By 
we denote an output signal which is calculated using the formula 

y^N = (16) 


Supposing that the matrix is given, to calculate the output signal it is 
necessary to perform N'^ complex multiplications and N{N — 1) complex ad¬ 
ditions. If the input signal is real, then the number of real multiplications will 
be equal to 27V^, and the number of real additions will be equal to 2N{N — 1). 

If we use the decomposition (14) of the matrix by calculating and 
will perform the matrix-vector multiplication from the right side to the left, 
the most time-consuming operations are multiplications of matrices and 
Vjv by the vector, because those matrices are not diagonal. If we interchange 
the columns of the matrix in the prescribed manner, we obtain a matrix 
Vat of special structure, which can be generated recursively. We will show it 
in Example It will allow to reduce the number of arithmetical operations 
by calculating the products of matrices and Vat by the vector. 
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Example 2 The matrices V^r for iV = 2,4 ,8 are as follows: 

1 -b -b 

-1-b] 
b 1 


V, = 


V4 = 


b -b^ -b 1 
b 1 -b -b^ 
b^ b 1 b 


Vs = 


1 -b b^ 

b - 6 ^ 63 
6 - 6 ^ -6 
62 -63 -62 


6 

62 

62 

63 


1 

6 

6 

62 


-6 

-62 

1 

6 


-6 

-62 

1 

6 

-62 


62 

-6 

-6 

1 

63 


-63 

62 

62 


62 

-6 

63 


-6 
1 

-62 
-6 -62 6 
62 -6 -&2 
-6 1 6 
-6 -62 -63 
1 6 62 


-63 -62 
6 -62 
62 6 

The matrix V 2 has some specihc structure. Now we consider the matrix V 4 . If 
in the matrix V 4 the second and fourth columns will be interchange and then 
the third and fourth columns will be interchange too, we obtain the following 
matrix: 

’ 1 -6 -6 62 ■ 

6 1 -62 -6 
6 -62 1 -6 
62 6 6 1 


V 4 = 


V 2 - 6 V 2 
6 V 2 V 2 


The matrix V 4 differs from the matrix V 4 only in the order of the columns. 
Therefore, the matrix V 4 can be obtained by post-multiplying the matrix V 4 
by the permutation matrix P 4 : 


V4 = V4P, 


where 


P 4 = 


Now we consider the matrix Vg 
of columns of this matrix: 

12 3 
1 8 4 


4i'4, 

'1 0 0 0 ‘ 

0 0 0 1 
0 10 0 
0 0 10 

If we perform the following permutation 


as a result we obtain the following matrix: 


Vg 


1 -6 -b 62 -6 62 62 -63 

6 1 -62 -6 -62 -6 63 62 

6 -62 1 -6 -62 63 -6 62 

62 6 6 1 -63 -62 -62 -6 

6 -62 -62 63 1 -6 -6 62 

62 6 -63 -62 6 1 -62 -6 

62 -63 6 -62 6 -62 1 -6 

63 62 62 6 62 6 6 1 


V 4 - 6 V 4 
6 V 4 V 4 
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As previously, we can write: 

V8=V8P8, 

where 

'1 0 0 0 0 0 0 0 ‘ 

00000001 
00010000 
00001000 
01000000 ■ 

00000010 
00100000 
_0 0 0 0 0 1 0 0 _ 

If we generalize the above considerations for = 2" we can write: 

(17) 

for TV = 2,4,, 2”. For N = 2 we can also write 

V2 = V2P2 = V2, 

where P 2 is an identity matrix of order two 

P2 = I2 

The permutation matrix Pjv of order iV = 2" can be obtained recursively 

from the permutation matrix Pjv /2 of order 2"“^ according to the following 

relation: 



(18) 


for TV = 4, 8 ,..., 2". Sn is the perfect shuffle permutation matrix of order 2", 
Jjv /2 is the counter-identity matrix of order N/2 and Ojv /2 is zero matrix. The 
perfect shuffle permutation is the permutation that splits the set consisting 
of an even number of elements into two piles and interleaves them. It can be 
written as follows: 

/l 2 3 4 ... 2n\ 


1 n 2 n + 1 ... 2n I ' 


For example 


Ss 


1 0 0 0 0 0 0 0 ' 
00001000 
01000000 
00000100 
00100000 ’ 
00000010 
00010000 
00000001 


'0 0 0 1 ' 
0 0 10 
0 10 0 
10 0 0 
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If we write the matrix Vjv as a product 'VjyPjsf the expression (14) will 
take the form: 


(19) 

The product PatA^P^ is a diagonal matrix, which has the same diagonal 
entries as the matrix but in different order and for a chosen parameter 
a it may be prepared in advance. If we denote this product multiplied by a 
factor 1/c" by A^: 


A% = —PmA%P 


NJ^N^N- 


the DFRHT algorithm (16) will take the following form: 


= VmA%v1^n 

where the matrix Vtv can be generated recursively: 


( 20 ) 


( 21 ) 


_ 

'1 -b' 

_ 

■ Vfc 

-bVk' 

V 2 = 

b 1 _ 

V2fc = 

bVk 

Vfc 


for fc = 2,4,...,2"-h 


( 22 ) 


3 Taking advantages of the particular structure of the matrix V 


N 


The most time-consuming operations by calculating the DFRHT transform 




according to (21) are multiplications of matrices 'Vjy and Vjv by the vector. 
Since in the matrix Vjv occur only following powers of 6 : 5", 6 "“^,..., 6° = 1 
we can write this matrix as follows: 


(0) 


V N = Ajy 


- 6A^^ + 5"a1"^ 


‘■iv 


■b^A 


(n) 

N 


(23) 

In the Figure it was shown the way of calculating the matrix-vector product 
ys = VgXg, using the expression (23). In this paper, the graph-structural 


models and data flow diagrams are oriented from left to right. Straight lines in 
the figures denote the operation of data transfer. We use the usual lines without 
arrows specihcally so as not to clutter the picture. Note that the circles in this 
figure shows the operations of multiplication by a number inscribed inside a 
circle. In turn, the rectangles indicate the matrix-vector multiplications by 
matrices 

Although it may seem strange, we will see that such an operation allows to 
reduce the number of multiplication and additions by multiplying the matrix 
Vat by a vector. It should be noted that because of the recursive relation (22) 


between matrices Vat and VAr/ 2 , the following recursive relation between the 


matrices A 


(k) A 
■N > ^N/2 


and aI^ occurs: 


'-N/2 


a(o) - 


,( 0 ) 


^N/2 ^N/2 
N/2 ^n/2 


0 


= I 


N 
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for fc = l, 2 , — 1 and n=log 2 iV, where 

= I2, 


a(o) - 
-^2 — 


1 0 
0 1 


^2 — 


0 -1 

1 0 


To clarify our idea we show the explicit form of expressions (23) and (24) for 
A = 2, A = 4 and A = 8 in an Example 


Example 3 

V 2 = + 6 A^^^ 

where the matrices A^^^ and A^^^ are presented above. 


where 


where 


V4 = A® + 6Ai^’ + b^A, 


( 1 ) 

*-4 


( 2 ) 

4 


a(o) _ 

A 4 — 


A 4 — 


A( 2 ) _ 
A 4 — 


10 0 0 
0 10 0 
0 0 10 
0 0 0 1 

0 - 1-1 0 
10 0-1 

10 0-1 

0 110 

0 0 0 1 

0 0-10 
0-100 
10 0 0 


,( 0 ) 


02 

a(0) 


— 


a(i) _a(°) 

a(0) a(i) 


Vg = A^°^ + hA),"’ + h^A\^’ + h^A 


( 1 ) 


O 2 -A 
aW 


( 2 ) 


( 1 ) 




10000000 

01000000 

00100000 

00010000 

00001000 

00000100 

00000010 

00000001 


(3) 


.( 0 ) 


04 


04 A® 


- 


0 - 1 - 10-10 0 0 

10 0 - 10-10 0 

1 0 0 -1 0 0 -1 0 

0110000-1 
1 0 0 0 0 -1 -1 0 

0100100-1 
0010100-1 
0 0 0 1 0 1 1 0 


a(i) _a(°) 

■^4 -^4 

a(0) a(1) 
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-1 

-1 

0 

0 

0 

0 

1 

1 

0 

0 

-1 

0 

0 

-1 0 
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-1 

0 

-1 
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0 

- 

0 0 

0 

0 

0 

0 0 

-1 
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0 0 

0 

0 

0 

0 1 
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0 0 

0 

0 

0 

1 0 
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0 0 

0 

0 

-1 

0 0 
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0 0 

0 

1 

0 

0 0 

0 



0 0 

-1 0 

0 

0 0 

0 



0-10 

0 

0 

0 0 

0 



1 0 

0 

0 

0 

0 0 

0 



*(2) _ a(1) 

-^4 -^4 

a(i) 

■^4 -^4 


O 4 -Af) 
Af^ O4 


Now we will evaluate the number of arithmetical operations, which are nec¬ 
essary to calculate the matrix-vector product Yn = Vjyxjv. We note, that in 
a general case such an operation requires multiplications and N{N — 1) 
additions. Now we will calculate the numbers of m ulti plications and additions 
needed for this operation if we use the expression (23) for the matrix V n, i.e. 


_ A 


( 0 )_ 


Since the non-zero entries of matrices A^^, A^^, ■ • •, A^^ are only 1 and - 
1 , no multiplications are needed when calculating the matrix-vector products 

(k) 

A]y^'xjv. The only multiplications we have to perform are multiplications of the 
vectors A^^xjv by the powers of b: A^^xjv by b, A^^xjv by 6 ^,..., A^^xat 
by 6". Because the number b is constant and known, its powers b^, , 6" 

may be prepared in advance. Thus, the number of multiplication by calcu¬ 
lating the matrix-vector product Vjyxjv is equal to nN = NlogN. Let us 
examine the number of additions, we need to perform, when calculating the 
matrix-vector product 'Vn^n- The total number of additions consist of num- 

(k) 

ber of additions by calculating the matrix-vector products Aj^^xjv, and nN 


additions which are needed to calculate the sum of vectors: A^^x^, 6A^^xjv, 




JV ) • 


^"■A^^xtv. Since, according to (24), the matrices A^^ have spe- 


(k) 

cific structures, the products A]^ xjv can be obtained by subtracting the prod¬ 
ucts A^^^ x^"^^ 




'-N/2 ^Af/2’ (excluding products A^^Xjv and A^^x^v which can 


N/2 ■^N/2 twice smaller size and summing the products 


be obtained even in a simpler way). By x^y 2 = [* 0 : a^i, ■ ■ ■, Xf^/ 2 -iY’ we denote 

the first half of the input vector xjv and by x ^^2 = [xn/ 2 tXm/ 2 +i^ ■ ■ ■ > 

- the second half of this vector, as it was shown, for A = 8, in the Figure 
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for fc = 1,2 


It should be noted that the products Aj^j 2 ^ 7 v /2 '^ 7 V/ 2 ^!v /2 

to calculate both products A^^xjv and A^'''^^X 7 v. For example the products 
A 4 °^X 4 ^^ and are used to calculate Ag^^xg and Ag^^xg. It allows to 

reduce the number of additions, because the some products are used twice. 

(fc) (1) 

Of course, this procedure can be repeated and each of products 

Aj^^^^x^y 2 J -^!v/2 ^7V/2 Calculated by summing (subtract¬ 

ing) products of twice smaller size and so on. It can be continued until cal¬ 
culating products of matrices and A^^^ by two-element sub-vectors of 
the vector x^r. The whole process of going down by calculating the product 
Yn = Vatxtv is presented, for A = 8, in the Figure]^ 


The expression (25) can be also written as the matrix-vector product, as 


follows: 


YN = VatXat 


CNx{n+l)N'^{n+l)N-^{n+l)NxN^N 


(26) 


where 
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CNx(n+l)N — llx(n+l) ® liV• 

The symbol 0 denotes the Kronecker product of matrices, and lix(n+i) is the 
matrix (row vector) whose all entries are equal to 1. The matrix A(„_|_i)jvxAr 
is responsible for multiplications of the matrices A^^ by the 

input vector xjv, the matrix - for multiplications of those products 

by the proper powers of b, and the matrix CnxuN - for aggregation of results. 
The process of going down by calculating the product V n^n, which has been 
presented in Figures and can be also described in the therm of matrices 
product. It means factorisation of the matrix A.(n+i)NxN into the product of 
n matrices 


A(„+i)ArxAf — A(„_|_i)^xnJvA„7Vx(n-l)Ar • ■ ■ A^ 


2NxN 


(27) 


The matrices which occur on the right side of expression (271 have the following 
forms: 

-^( 0 ) 

A-2NxN = lAr/2 ® (28) 

[-^ 2 x 2 


where 


A2X2 = Af 0 [1] 0 Ii = Af, A^',^2 = A^'^ 0 [1] 0 Ii = A^^h 


.( 1 ) 


^3Nx2N 


= 1 


N/4 ' 


^( 0 ) 

-^4x8 

-^4x8 + -^4x8 

-^4x8 


(29) 


where 

A ®8 = a(°^ ® [1 0 ] 0 I 2 , = A^^^ 0 [0 1 ] 0 I 2 

and the matrix denotes the matrix A^^^^g which columns were circu¬ 
larly shifted by 2 positions to the right, and the matrix denotes the 

matrix A^^^^g which columns were circularly shifted by 2 positions to the left. 
The last matrix A.(^n+i)NxnN is defined as 


A(„+i)ArxnAf — ItV/TV 


X( 0 ) 

■^N xnN 

^NxnN -^NxnN 


—{0)an-l)N/2^) (l)(N/2^) 

^NxnN + ^NxnN 


^NxnN 


(30) 


where 

AATxnAr = A 2 ^ 0 [1 0 ... 0](8)Ijv/2j = A 2 ^ ® [0 0 ... l]®Ijv/ 2 . 
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Using the expression (27) the algorithm (26) of calculating the product V nxn 


can be written as follows: 
Yn = C 


Nx(n+l)N'B(n+l)N-^(n+l)NxnN-^nNx(n-l)N ■ ■ ■ ^2NxN^N (31) 


The expression (31) allows for evaluating the total number of additions, which 


are needed to calculate the matrix-vector product V^tXat. We assume that 
the input vector xjv is real-valued. Each of matrices A(fe_|_i)AfxfciVj for k = 
1,2,... ,n, is the direct sum of iV/2^ identical blocks and the single block is 
the vertical concatenation ol k + 1 matrices. The firs, indicated by A 2 fcxfc 2 *=i 

and the last, indicated by A. 2 '‘xk 2 >‘^ do not need any additions or subtrac¬ 
tions by multiplying them by a vector. The k — 1 others matrices, which are 

sums of A 2 fexfc 2 '= ^'od A 2 fexfc 2 '=: after circularly shifting their columns, so mul¬ 
tiplying each of them by a vector needs 2^ additions. To calculate the prod¬ 
uct A(„+i)Arx„ArA„Arx(„-i)Ar ■.. A 2 Nxn^n we have to perform J2k=i |v(fc- 
1)2^ = Nn{n — l)/2 additions. The product of the matrix B(„_|_i)Ar by a vec¬ 
tor do not need any additions and the product of the matrix CArx(n-i-i)Af by a 
vector needs nN additions. Thus t he to tal number of additions by calculating 
the products VjvXat, according to (31), is equal to Nn{n + l)/2. 

Example shows the explicit form of the algorithm (31) with all occurring 
in it matrices for N = 8. 


Example 4 The algorithm (31) of calculating the product of matrix Vg by the 
vector Xg is as follows: 


ys — VgXg — Cgx32B32A32x24A24xl6Ai6x8Xg, 


where 


^16x8 — i4 ' 


^( 0 ) 

-^ 2 x 2 

-^2x2 


^( 0 )( 2 . 

^4x8 


1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 
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0 6 0 0 
0 0 62 0 
0 0 0 63 



® Is , 








^8x32 — llx4 ® Is - 

It is easy to check that in this case the total number of addition is equal to 
48 and the number of multiplications is equal to 24 (we can see it also in the 
figure [^. 


4 The novel DFRHT algorithm 


Now we return to the DFRHT algorithm (21). According to (23) the matrix 


VAT can be written as the sum of the matrices , 


, A^^ with CO- 
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efficients 1, The transposed matrix can be written as the sum 

of the transposed matrices , A^^^ with the same coefficients 


1 , 6 ,..., 6 ": 


_ A (0)^ 
^ N — 


- 6 aW^ + 62 . 


6 ”A 


(n)T 

N 


(32) 


( 0 ) A ( 2 ) 


N > 


'■N ’ ■ ■ 


are symmetric and 


Since the matrices with the even indexes A 
the matrices with the odd indexes are asymmetric the expression (32) can be 
written as follows: 


= A)^' - 6A« + ... + (-1)”6”A 


(n) 

N 


(33) 


According to (261 the matrix V n from expression (231 can be transformed 
into the product 


Vat — CArx(n+l)AfB(„_|_i)ArA(„_|_i)ArxAr 


(34) 




SO the matrix V^ may by also transformed from (331 into the following prod¬ 
uct: 


V^ = C, 


Afx(n-|-l)ArB(„_|_i)ArA(„_|_i)Arx7V 


(35) 


where the matrix CArx(Ti+i)Ar is defined as follows: 


CiVx(n-|-l)Ar — llx(n-l-l) ® In 


and 


Llx(n-l-l) 


= [1 -1 ... (- 1 )"]. 


The others matrices in the expression (35) are the same as that in the expres¬ 


sion (26). Taking into account the decompositions (34) and (35) of matrices 




VAT and respectively the DFRHT algorithm (21) will take the form: 


— CNx(n+l)N'^{n+l)N-^{n+l)NxN^N^Nx(n+l)N'^(n+l)N-^{n+l)NxN^N 

(36) 


where the matrix A(„+i)ArxJV is decomposed according to (27). For example, 
for A = 8 this algorithm will take the following form: 


yt'’ = C 8 X 32 B 32 A 32x24A24xl6Ai6x8AgC8x32B32A32x24A24xl6Ai6x8X8. 


Figure]^ shows a data flow diagram of the algorithm for 8 point DFRHT. 
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Fig. 4 


Data flow diagram of the DFRHT algorithm 



for iV = 8 


5 Discussion of computational complexity 

Let us assess the computational complexity in term of numbers of multiplica¬ 
tions and additions required for DFRHT calculation. Calculation of the dis¬ 
crete fractional Hadamard transform for a real-valued vector xjv of length 
N = 2'^, assuming that the matrix defined by Q is given, requires = 
2^" multiplications of a complex number by a real number and N{N — 1) = 
2’^(2" — 1) complex additions. Each multiplication of a complex number by a 
real number needs two real multiplications and each addition of two complex 
numbers requires two real additions. Hence the numbers of real multiplications 
and real additions required for computing the DFRHT using the naive method 
are equal to 2^"’+^ and 2"’+^ (2" — 1) respectively. 

Let us now evaluate the computational complexity of the DFRHT realiza¬ 
tion with the help of the procedure ( |36| ). As it was discussed in the section 
if we use the factorized representation of the matrices and V^y, calculat- 
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—T 

ing the product of the real-valued matrix and the real-valued vector xjv 
requires nN real multiplications and Nn{n + l)/2 real additions. As a result, 
we again obtain the real-valued vector. Then there is computed the product of 
the complex-valued diagonal matrix and the real-valued vector calculated 

previously (we assume that for a predetermined parameter a, the diagonal 
elements of this matrix were calculated in advance). The calculation of this 
product requires 2N real multiplications. The resulting complex-valued vector 
is then multiplied by the factorized matrix V^v- This operation requires 2Nn 
real multiplications and Nn(n -I- 1) real additions. The total numbers of arith¬ 
metic operations to compute DFRHT of size 2" using our new algorithm are 
N{3n + 2) real multiplications and 3Nn{n + l)/2 real additions. It is easy to 
check that even for small n the numbers of arithmetic operations required for 
realization of the proposed algorithm are several times less than in the naive 
method of computing. 

Tables and display the numbers of multiplications and additions re¬ 
quired for the DFRHT transform implementation of the real-valued input 
signal of the length = 2”. These numbers were calculated for three methods 
of the transform implementation: the direct multiplication of the DFRHT ma¬ 
trix by a vector of the input data, calculation according to authors’ algorithm 
described in the work El, and according to the algorithm ( [M| ) proposed in 
this article. It is easy to check that for n > 5 the number of arithmetic oper¬ 
ations, required for DFRHT transform realization according to the proposed 
algorithm, is smaller than in the other two methods of DFRHT computing. 


Table 1 Numbers of multiplications for mentioned algorithms 


Af = 2" 

direct method 

method [Tl] 

proposed algorithm 

2 

8 

6 

10 

4 

32 

18 

32 

8 

128 

54 

88 

16 

512 

162 

224 

32 

2048 

486 

544 

64 

8192 

1458 

1280 

128 

32768 

4374 

2944 

256 

131072 

13122 

6656 

512 

524288 

39366 

14848 

1024 

2097152 

118098 

32768 


6 Summary 

The article presents the novel algorithm for the DFRHT performing. The 
algorithm has a much lower computational complexity than the direct way 
of the DFRHT implementation. The computational procedure for DFRHT 
calculating is described in Kronecker product notation. The Kronecker prod¬ 
uct algebra is a very compact and simple mathematical formalism suitable 
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Table 2 Numbers of additions for mentioned algorithms 


Af = 2” 

direct method 

method 

proposed algorithm 

2 

4 

5 

6 

4 

24 

25 

36 

8 

112 

95 

144 

16 

480 

325 

480 

32 

1984 

1055 

480 

64 

8064 

3325 

1440 

128 

32512 

10295 

4032 

256 

130560 

31525 

10752 

512 

523264 

95855 

69120 

1024 

2095104 

290125 

168960 


for parallel realization. This notation enables us to represent adequately the 
space-time structures of an implemented computational process and directly 
maps these structures into the hardware realization space. For simplicity, we 
considered the synthesis of a fast algorithm for the DFRHT calculation for 
N = 2^. However it is clear that the proposed procedure was developed for 
the arbitrary case when the order of the matrix is a power of two. 
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